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ABSTRACT 

Singularly  perturbed  models  of  the  form  x(t )  =  A(e)x(t),  with  A(e)  analytic  at  0,  non¬ 
singular  for  e  €  (0,co]  and  singular  at  e  =  0,  arise  naturally  in  various  problems  of  systems  and 
control  theory.  Under  a  so-called  multiple  semi-simple  null  structure  or  MSSNS  condition  on 
A(e),  the  eigenstructure  of  this  matrix  has  a  multiple  scale  property  that  allows  the  asymptotic 
eigenstructure  of  the  matrix  to  be  studied  via  reduced-order  matrices  associated  with  the  sep¬ 
arate  scales.  Under  a  stronger  multiple  semi- stability  or  MSST  condition,  this  eigenstructure 
decomposition  translates  into  a  time-scale  decomposition  of  the  solution  x(t)  of  the  system. 
This  paper  is  aimed  at  illuminating  the  MSSNS  and  MSST  conditions.  Using  ideas  from  an  al¬ 
gebraic  approach  that  we  have  developed  for  the  study  of  singularly  perturbed  systems,  we  show 
(among  other  results):  that  the  Smith  decomposition  of  A(e)  permits  transformation  to  a  form 
in  which  the  scales  become  explicit  and  the  computations  become  transparent  and  simple;  that 
this  form  allows  us  to  identify  perturbations  of  A(e)  that  all  have  the  same  scale-based  decom¬ 
positions;  that  the  eigenstructure  of  A(e)  does  indeed  display  multiple  scales  under  MSSNS, 
and  that  the  eigenvalues  in  this  case  can  be  approximated  using  reduced-order  calculations; 
that  A(e)  has  MSSNS  if  and  only  if  the  orders  of  its  invariant  factors  equal  the  orders  of  its 
eigenvalues;  that  this  happens  if  and  only  if  the  orders  of  its  invariant  factors  and  principal 
minors  are  related  in  a  specified  way;  and  that  A(e)  has  MSST  if  and  only  if  it  is  Hurwitz  for 
e  €  (0,eo]  and  has  MSSNS  and  satisfies  the  requirement  that  the  order  in  e  of  the  real  part  of 
every  eigenvalue  is  not  greater  than  the  order  of  its  imaginary  part. 
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1.  INTRODUCTION 


Various  problems  in  systems  and  control  theory  give  rise  to  singularly  perturbed ,  linear, 
time  invariant,  n-th  order  models  of  the  form 

i(i)  =  A(e)z(f)  (1) 

with  the  matrix  A(e)  being  n  x  n,  analytic  at  0,  nonsingular  1  for  e  e  (0,co]  and  singular  at 
e  =  0;  see  the  recent  books  [l],  [2]  and  survey  [3],  Such  models  are  found  in,  for  example,  studies 
of  root  locus  behavior  under  high  gain  feedback  [4]-[6] ,  in  the  analysis  of  power  systems  [7],  and 
in  the  context  of  Markov  models  [8],  [9],  The  small  positive  parameter  e  in  these  respective 
contexts  may,  for  example,  represent  the  inverse  of  a  high  gain  in  a  feedback  controller,  the 
inverse  of  the  large  rotational  inertia  of  some  generator,  or  a  small  failure  rate  or  inspection 
rate  in  a  reliability  model. 

A  key  objective  in  the  study  of  singularly  perturbed  models  is  to  express  the  behavior  of 
(1)  as  a  perturbation  of  its  behavior  for  e  =  0.  In  particular,  since  A(e)  becomes  singular  at 
e  =  0,  attention  focuses  on  behavior  associated  with  those  eigenvalues  of  A(e)  that  tend  to  0  as 
e  tends  to  0,  i.e.  on  the  “zero  group”  of  eigenvalues  of  A(e),  [10].  A  prototype  for  results  on 
singularly  perturbed  systems  is  provided  by  the  following  example. 

Example  1  Suppose 


where  the  A<y  are  constant  matrices.  If  An  is  nonsingular,  then  it  is  not  hard  to  show  that 
the  eigenvalues  of  A(e)  fall  into  two  separate  groups:  a  fast  group  that  is  asymptotically  domi¬ 
nated  by  the  eigenvalues  of  An,  and  a  slow  group  —  the  zero  group  —  that  is  asymptotically 
dominated  by  the  eigenvalues  of  eA22 ,  where  A22  is  the  Schur  complement  of  An  in  the  matrix 


i.e., 

A22  =  A22  —  A2i  A111  Ai2  (4) 

Note  that  our  standing  assumption  on  the  nonsingularity  of  A(e)  away  from  0  guarantees  that 
A  in  (3)  is  nonsingular.  With  this  assumption,  it  is  easy  to  see  that  the  assumed  nonsingularity 
of  An  implies  the  nonsingularity  of  A22. 

The  role  of  An  in  approximating  the  fast  group  of  A(e)  is  easy  to  understand  on  the 
basis  of  simple  continuity  arguments,  since  the  nonzero  eigenvalues  of  A(0)  are  those  of  An. 

1  Our  results  can  be  extended  quite  easily,  see  [14],  [15],  to  the  case  where  A(e)  has  semi-simple 
null  structure,  i.e.  where  the  eigenvalue  at  0  has  geometric  multiplicity  equal  to  its  algebraic  multiplicity, 
but  we  avoid  this  case  for  simplicity  here. 
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The  appearance  of  A22  requires  a  little  more  motivation.  For  this,  apply  to  A(e)  a  similarity 
transformation  defined  by  the  matrix 

rW=(o  A“/‘12)  (s) 

thereby  obtaining 

r(f)AWr-‘W=(A'^"  e£j  M 

where  the  *’s  denote  constant  matrices.  Since  the  coupling  between  the  nonsingular  diagonal 
blocks  that  is  induced  by  the  off-diagonal  blocks  appears  smaller  in  (6)  than  in  (1),  one  might 
now  be  more  inclined  to  believe  the  above  claim  regarding  approximation  of  the  fast  and  slow 
groups  of  eigenvalues  of  A(e).  A  rigorous  proof  can  be  obtained  by  applying  well  established 
arguments  to  the  matrix  in  (6),  for  example  a  block  version  of  Gerschgorin’s  theorem,  [11],  or 
the  arguments  in  [12],  which  we  follow  in  Section  3. 

Further  intuition  for  this  separation  result  will  be  obtained  when  we  continue  the  example 
later  in  this  section.  The  conclusion  for  now  is  that  knowing  the  eigenvalues  of  An  and  A2s 
allows  one  to  determine  the  dominant  behavior  of  the  eigenvalues  of  A(e),  and  hence  to  obtain 
asymptotically  good  approximations  of  them  via  reduced-order  calculations.  It  should  already 
be  evident  from  the  sketch  above  that  the  separation  results  are  unchanged  if  the  A*y  in  (2) 
have  higher  order  terms  appended  to  them.  This  fact  will  be  elaborated  on  in  Section  2.  •• 

The  special  form  of  A(e)  in  (2)  made  the  computations  above  quite  transparent,  while  the 
condition  that  An  be  invertible  made  the  computations  possible.  The  invertibility  condition 
is  equivalent  to  this  special  A(e)  having  what  is  termed  multiple  semi-simple  null  structure  or 
MSSNS,  [8],  [13].  The  MSSNS  condition  for  a  general  A(e)  plays  a  fundamental  role  in  studies 
of  its  asymptotic  eigenstructure,  and  a  primary  objective  of  the  present  paper  is  to  illuminate 
this  condition  by  presenting  certain  equivalent  forms  of  this  condition  and  by  examining  certain 
consequences  of  MSSNS.  We  begin  in  Section  2  with  a  review  of  results  in  [14],  [15]  that  show 
how  the  Smith  decomposition  of  a  general  A(e)  yields  a  similarity  transformation  that  brings 
this  matrix  to  a  special  form  —  which  we  term  the  (reduced)  explicit  form  —  in  which  MSSNS 
is  easily  checked  and  scale-based  decomposition  easily  carried  out.  This  explicit  form  is  an 
extension  of  the  one  in  Example  1  above,  and  checking  it  for  MSSNS  reduces  to  checking 
invertibility  of  a  sequence  of  Schur  complements.  The  explicit  form  also  allows  us  to  identify 
perturbations  of  A(e)  that  have  the  same  scale-based  decompositions. 

Building  on  the  explicit  form,  Section  3  deduces  certain  consequences  of  MSSNS  and 
certain  equivalent  tests  for  it.  In  particular,  it  is  shown  that  the  eigenstructure  of  A(e)  displays 
multiple  scales  under  MSSNS,  and  that  the  eigenvalues  in  this  case  can  be  approximated  using 
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reduced-order  calculations.  It  is  also  demonstrated  that  A(e )  has  MSSNS  if  and  only  if  the 
orders  of  its  invariant  factors  equal  the  orders  of  its  eigenvalues. 

Since  eigenvalues  are  determined  by  principal  minors  while  invariant  factors  are  deter¬ 
mined  by  all  the  minors,  it  is  evident  that  MSSNS  must  imply  some  relation  between  these 
two  sets  of  minors.  We  show  in  Section  4  that  in  fact  the  MSSNS  condition  is  satisfied  if  and 
only  if  the  orders  of  the  invariant  factors  and  principal  minors  of  A(e)  are  related  in  a  specified 
way.  The  insights  provided  by  these  results  can  be  applied,  as  outlined  in  that  section,  to  the 
problem  of  scaling  a  matrix  A(e)  via  a  diagonal  similarity  transformation  to  induce  MSSNS  in 
it. 

Many  of  the  results  on  MSSNS  in  Sections  3  and  4  turn  out  to  echo  results  in  [5].  The 
development  in  [5]  is  in  the  context  of  asymptotic  root  loci  of  systems  under  high-gain  feedback, 
and  in  principle  a  mapping  can  be  made  between  the  formulation  there  and  ours  here.  Though 
our  route  to  results  of  interest  for  the  system  (1)  is  more  direct,  the  treatment  in  [5]  is  of  value 
in  rounding  out  our  treatment,  and  the  interested  reader  is  encouraged  to  examine  that  paper 
as  well. 

The  illustration  in  Example  1  involved  eigenstructure  decomposition.  One  can,  with 
further  assumptions,  go  beyond  this  to  a  time-scale  decomposition  of  the  solution  x(t)  of  (l), 
as  shown  next. 

Example  1,  Continued  Consider  the  system  (1)  with  A(e)  defined  as  in  (2),  so  that 

^12  \  ( xi  (t)\  /7\ 

\x2{t)J  Ve^2i  cA22J\x2(t)J  '  ' 

Much  of  the  systems  and  control  literature,  see  [l]-[3],  focuses  on  this  prototype  example  (or 
one  that  can  be  obtained  from  it  by  simply  changing  the  independent  variable  from  t  to  r  —  et) . 
Take  An  to  be  nonsingular,  as  before.  If  e  —  0,  then  x2(t)  evidently  remains  constant  at  its 
initial  value  x2(0),  while  a?i(t)  is  governed  by  the  fast  group  of  eigenvalues  of  A(0),  namely  the 
eigenvalues  of  An-  Under  the  assumption  that  An  is  Hurwitz  and  not  just  nonsingular,  i.e. 
that  all  its  eigenvalues  lie  in  the  open  left  half  of  the  complex  plane,  xx  (t)  settles  down  after 
a  fast  transient  to  a  steady  state  value  that  is  readily  calculated  by  setting  xL  (t)  =  0  in  (7): 
xi(oo)  =  -A~1A12X2(o o)  =  -A^1  A12x2(0). 

Under  the  same  assumption  of  Hurwitz  An,  the  behavior  for  small  nonzero  e  is  only 
slightly  different  from  this:  x2(t)  now  varies  slowly  instead  of  being  constant,  so  that  the 
behavior  of  xx (t),  after  a  fast  transient  governed  approximately  by  the  eigenvalues  of  An,  is 
well  approximated  by  —  A^1  Ai2x2(t).  Using  this  latter  approximation  for  Xi(t)  in  the  equation 
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for  x2  (t)  yields  the  following  approximate  governing  equation  for  x2  ( t ) : 


x2{t)  «  e(A22  —  A2iA1^ Ai2)x2(t)  —  eA22x2 (t) 


(8) 


The  appearance  of  the  matrix  eA22  here  is  consistent  with  what  was  claimed  earlier,  namely 
that  the  slow  group  of  eigenvalues  of  A(e)  is  well  approximated  by  the  eigenvalues  of  eA22.  If 
we  assume  that  A22  is  also  Hurwitz,  (8)  yields  a  uniformly  good  approximation  to  the  solution 
over  an  infinite  interval.  To  summarize  all  this  more  precisely,  what  can  be  shown  under  the 
above  assumptions  is  that 


*i(0 

x2(t) 


xlf  (t)  +  xu(et)  +  0(€) 
x2,(et )  +0(e) 


t  >  0 


where 


Xif{t)  =  AuXlf{t),  xlf{0)  =  Xi(0)  -  xx.(O) 
xi,(t)  =  -A~11A12x2s  ( t ) 
x2,  ( t )  =  A22  x2,  ( t ) ,  x2,  (0)  =  x2  (0) 


(9a) 


(96) 

(9c) 

m 


and  the  subscripts  /  and  s  respectively  denote  “fast”  and  “slow”  parts  of  the  solutions. 


The  Hurwitz  conditions  above  are  equivalent  to  having  the  special  A(e)  in  the  example 
satisfy  a  so-called  multiple  semi- stability  or  MSST  condition,  [8],  [13] .  The  MSST  condition  for 
a  general  A(e)  plays  a  fundamental  role  in  the  time-scale  decomposition  of  singularly  perturbed 
systems,  see  [8]  and  [13]-[15],  The  latter  reference,  [15],  contains  particularly  strong  results  on 
the  connection  between  MSST  and  time-scale  decompositions.  In  the  present  paper,  our  focus 
is  not  on  the  time-scale  decomposition  itself  but  on  features  of  the  MSST  condition  and  on  its 
relation  to  other  system  properties.  The  MSST  condition  is  introduced  along  with  MSSNS  in 
Section  2,  where  the  test  for  MSST  for  systems  in  explicit  form  is  presented.  We  then  show  in 
Section  3  that  A(c)  has  MSST  if  and  only  if  it  is  Hurwitz  for  e  €  (0,eo]  and  has  MSSNS  and 
satisfies  the  requirement  that  the  order  of  the  real  part  of  every  eigenvalue  is  not  greater  than 
the  order  of  its  imaginary  part.  This  result  suggests  the  role  of  MSST  in  obtaining  time-scale 
decompositions . 
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2.  DECOMPOSITION  CONDITIONS  FOR  A(e)  IN  EXPLICIT  FORM 
2.1  Definitions  of  MSSNS  and  MSST 

Before  presenting  our  algebraic  approach  to  scale-based  decompositions  and  the  condi¬ 
tions  that  enable  these  decompositions,  we  shall  review  the  definitions  of  MSSNS  and  MSST 
presented  in  [8],  [13].  First  note  that  a  matrix  is  said  to  have  semi-simple  null  structure  or 
SSNS  if  the  geometric  multiplicity  of  the  eigenvalue  at  0  (when  it  is  present)  equals  its  alge¬ 
braic  multiplicity;  it  is  said  to  be  semi-stable  if  it  has  SSNS  and  all  its  nonzero  eigenvalues  have 
negative  real  parts. 

Since  A(e)  is  analytic  at  0,  it  has  a  Taylor  expansion: 

oo 

A(e)=X>PF\p  (10) 

P=  0 

Let  M^e)  be  the  total  projection  for  the  zero  group  of  eigenvalues  of  A(e),  i.e.  the  projection 
onto  the  corresponding  eigenspace,  [10].  If  Fl0  in  (10)  has  SSNS,  then  it  can  be  shown  that 

A2(e)  =  M1(e)A(e)/c  (11a) 

also  has  a  Taylor  expansion,  of  the  form 

OO 

A2(e)  =  ]>>F2p  (116) 

p=0 

If  F2 o  also  has  SSNS,  we  can  similarly  define 

oo 

A3(e)  =  M2(e)A2(e)/e=  ^  ep  F3p  (12) 

p= o 

where  M2(e)  is  the  total  projection  for  the  zero  group  of  eigenvalues  of  A2(e).  This  process  can 
be  continued,  terminating  when  Fm0  does  not  have  SSNS  or  when 

p{Fio)  +  •  •  -p(Fmo)  —  n  (13) 

where  p(F)  denotes  the  rank  of  F  and  n  is  the  order  of  the  system  (1),  i.e.  the  dimension  of 

A(e). 

Now  A(e)  is,  by  definition,  said  to  satisfy  MSSNS  if  the  above  construction  can  proceed 
all  the  way  to  the  stage  m  at  which  (13)  is  satisfied,  with  Fi0 Fm0  all  having  SSNS. 
If  in  addition  F10,  ■  ■  ■ ,  Fm0  are  all  semi-stable,  then  A(e)  is,  by  definition,  said  to  satisfy 
MSST. 

The  significance  of  these  two  conditions  emerges  from  the  results  in  [8],  [13]-[15].  In 
particular,  [15]  shows  that  (1)  has  what  is  termed  a  strong  time-scale  decomposition  if  and  only 


if  it  satisfies  MSST.  If  the  system  satisfies  only  MSSNS  but  is  stable  for  positive  e,  then  it  may 
still  be  possible  to  obtain  what  [14]  and  [15]  term  am  extended  time-scale  decomposition.  In 
any  case,  it  is  always  possible  under  MSSNS  to  obtain  an  eigenstructure  decomposition,  as  will 
be  shown  in  Section  3. 

2.2  Transformation  to  Explicit  Form  Using  the  Smith  Decomposition  of  A(e) 

The  starting  point  for  our  results  is  the  Smith  decomposition  of  A(c).  Note  that  the  entries 
of  A(e)  are  elements  of  the  (local)  ring  of  functions  analytic  at  e  =  0,  i.e.  functions  expressible 
as  Taylor  series  in  e.  This  ring,  which  we  shall  denote  by  W,  is  a  Euclidean  domain  [16],  with 
the  degree  of  an  element  being  the  order  of  the  first  nonzero  term  in  its  Taylor  expansion  — 
the  degree  of  a,e*  +  a,+  1e,  +  1  +  . .. ,  yt  0  is  thus  i.  The  units  in  this  ring  are  precisely  those 
elements  that  have  degree  0,  i.e.  those  that  are  nonzero  at  e  =  0.  A  square  matrix  17(e)  over 
this  ring,  i.e.  one  with  entries  from  this  ring,  is  termed  unimodular  if  det  U (e)  is  a  unit,  or 
equivalently  if  det  17(0)  ^  0,  or  equivalently  if  U~1(e)  is  also  over  this  ring. 

It  now  follows  from  well  known  results  on  Smith  forms  of  matrices  over  Euclidean  do¬ 
mains,  [16],  that  the  n  x  n  matrix  A(e)  has  the  decomposition 

A(e)  =  P(e)D(e)Q(e)  (14a) 

where  P(e)  and  Q(e)  are  unimodular  matrices  over  W  and 

D(e)  =  block  diagonal  ( Ikl ,  elki ,...,  em~1Ikm  )  (146) 

with  Ik.  denoting  an  identity  matrix  of  dimension  k}  x  fcy,  with  k}-  =  0  corresponding  to  absence 
of  the  y-th  block,  and  with  km  ^  0.  The  indices  k3,  and  hence  D(e),  are  unique,  though  P(e) 
and  Q(e)  are  not. 

The  above  indices  actually  capture  the  invariant  factor  structure  of  A(e);  in  the  pole-zero 
language  of  [17],  which  is  devoted  to  rational  matrix  structure  but  has  much  that  applies  to 
local  rings,  the  indices  represent  the  zero  structure  of  A(e)  at  e  =  0.  The  t'-th  invariant  factor, 
by  definition,  is  the  ratio  of  the  greatest  common  divisor  (gcd)  of  all  i  x  i  minors  of  A(e)  divided 
by  the  gcd  of  all  (t  —  1)  x  (t  —  1)  minors,  with  the  first  invariant  factor  being  defined  as  the  gcd 
of  all  the  lxl  entries.  Since  the  gcd  is  only  defined  up  to  a  unit,  we  can  always  represent  it  in 
the  ring  W  by  an  element  of  the  form  er  for  some  nonnegative  integer  r,  so  the  invariant  factors 
are  also  of  this  form.  Now  what  the  Smith  decomposition  in  (14)  captures  is  the  fact  that  A(e) 
has  k}-  invariant  factors  of  the  form  cy_1,  i.e.  k3  invariant  factors  of  (degree  or)  order  j  —  1. 

The  decomposition  (14)  allows  us  to  similarity  transform  the  given  model  (1)  to  an 
equivalent  model  in  which  the  potential  for  scale-based  decomposition  is  considerably  more 
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transparent,  and  where  the  analysis  is  as  direct  as  in  Example  1.  For  this,  define 

y(t)  =  P~1(c)x(t)  (15) 

so  that,  combining  (1)  and  (14),  one  gets  the  explicit  form  model 

y{t)  -  D(e)A(e)y{t)  (16a) 

or 


(I* 

0 

0  A 

Mu(e) 

■4^X2  (e) 

...  Alm  (e)  y 

2/2(4 

— 

0 

0 

A21  (e) 

A22  (f) 

■  ■  -  A2m  (e) 

y2(4 

U m(t)J 

V  o 

0 

■rffl  —  1  T  . 

■■  e  lkm  J 

i(e] 

1  4^2  (e) 

. . .  Amm  (e)  J 

<2 /m(0> 

(166) 

where  the  matrix  A(e)  is  unimodular  and  given  by 

A(e)  =  Q(e)P(e)  (16c) 

(We  have  dispensed  with  overbars  for  the  entries  of  A(e )  in  (16b)  so  as  to  simplify  the  notation.) 
Since  the  unimodular  transformation  matrix  P(e)  in  (15)  is  finite  and  nonsingular  at  e  =  0,  the 
asymptotic  behavior  of  (1)  may  be  retrieved  without  difficulty  from  that  of  the  equivalent  system 
(16).  It  is  also  easy  to  see  from  Section  2.1  that  such  a  unimodular  similarity  transformation 
does  not  alter  MSSNS/MSST.  (For  more  on  such  “analytic  similarity”  of  matrices,  see  [18].) 

We  can  go  still  further.  Examination  of  (I6a,b)  might  suggest  that  the  potential  for 
scale-based  decomposition  is,  at  least  under  appropriate  conditions,  completely  displayed  by 
D(e),  with  the  unimodular  matrix  A(e)  not  contributing  anything.  This  thought  would  lead 
one  to  examine  a  model  in  which  A(e)  is  replaced  by  A(0): 

’z{t)  =  D{e)A{0)z(t)  (17) 

We  shall  term  this  the  reduced  explicit  form  of  (1).  A  key  result  in  [14]  (see  also  [15]), 
obtained  by  carrying  out  the  computations  of  Section  2.1  in  detail  for  A(e)  given  by  (14), 
is  that  the  models  in  (1),  (16)  satisfy  MSSNS  (respectively  MSST)  if  and  only  if  the  above 
reduced  explicit  model  satisfies  MSSNS  (respectively  MSST).  Furthermore,  the  eigenstructure 
decomposition  of  (17)  obtained  under  the  MSSNS  condition  is  an  eigenstructure  decomposition 
of  (16)  and  (1);  this  will  be  discussed  in  Section  3.  Similarly,  if  MSST  is  satisfied,  then  a 
time-scale  decomposition  of  (17)  is  shown  in  [14],  [15]  to  be  a  time-scale  decomposition  of  (16) 
as  well,  and  to  yield  a  time-scale  decomposition  of  (l)  on  transformation  by  P(0). 

These  results  therefore  justify  our  discarding  all  e-dependent  terms  in  A(e),  if  our  interest 
is  only  in  testing  for  MSSNS/MSST  and  effecting  the  scale-based  decompositions  above.  This 
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is  a  major  simplification,  since  tests  of  the  MSSNS/MSST  conditions  in  (17),  and  computations 
associated  with  scale-based  decompositions  of  it,  are  as  transparent  and  direct  as  those  asso¬ 
ciated  with  the  model  in  Example  1  (which  is  already  in  reduced  explicit  form).  The  existing 
literature  on  systems  in  the  special  form  (7)  can  therefore  be  easily  applied,  with  straightforward 
extensions,  to  the  much  more  general  systems  described  by  (1),  once  the  Smith  decomposition 
(14)  of  A(e)  has  been  determined. 

Another  important  consequence  of  the  results  above  is  that  inferences  made  regarding  (1) 
on  the  basis  of  computations  with  (17)  also  hold  for  any  other  system  that  has  the  same  reduced 
explicit  form,  (17).  The  set  of  matrices  iA(f)  that  give  rise  to  the  same  reduced  explicit  form 
(17)  as  A(e)  is  precisely  given  by 

Ah.a(c)  =  R(€)P{€)D{e){A{0)  +  A{e)}p-1{e)R-l{e)  (18a) 

where  R(e)  is  any  unimodular  matrix  and  A(e)  is  any  matrix  over  W  such  that  A(0)  =  0.  This 
can  be  rewritten  as 

A*,  a  W  =  R(e)A(e)[Q-l(e)[Q(0)P(0)  +  A  (e))P~'(e)]R-'(e) 

=  R(e)A{e)UA{e)R~1{e)  (186) 

where  UA(c)  is  a  unimodular  matrix  such  that  UA(0)  =  I.  Each  A(e)  gives  rise  to  a  unique 
unimodular  matrix  satisfying  this  condition,  and  conversely.  The  set  of  interest  is  thus  precisely 
characterized  by  (18b). 

The  key  to  the  above  development  is  the  Smith  decomposition  of  A(e).  Computation  of 
the  Smith  form  is  nontrivial,  of  course,  and  involves  operations  comparable  to  those  required 
by  [8],  [13]  to  compute  the  Fip  in  (10)-(13)  (though  we  have  found,  for  examples  that  are  small 
enough  to  work  out  by  hand,  that  computation  of  the  Smith  form  is  decidedly  simpler  than 
the  necessary  operations  on  the  matrix  coefficients  of  Taylor  series).  However,  our  approach 
here  permits  the  analysis  of  (1)  to  be  separated  into  a  transformation  step,  which  involves 
determination  of  the  Smith  form  and  produces  the  reduced  explicit  form,  and  a  greatly  simpli¬ 
fied  decomposition  step.  In  contrast,  the  operations  of  transformation  and  decomposition  are 
interleaved  in  the  approach  of  [8]  and  [13], 

Partly  because  of  its  two-step  nature,  and  partly  because  of  the  algebraic  connections, 
our  approach  has  also  yielded  several  new  insights  and  produced  new  results,  for  example 
on  feedback  assignment  of  closed-loop  time-scales,  on  perturbations  that  preserve  time-scale 
decompositions,  on  e-dependent  amplitude  scaling  of  state  variables  to  induce  MSSNS,  and  so 
on,  see  [14],  [15]  and  results  in  this  paper. 
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2.3  Checking  MSSNS  and  MSST  in  the  Reduced  Explicit  Form 

To  simplify  notation  for  the  reduced  explicit  form,  denote  A(0)  of  (17)  simply  by  A,  and 
its  entries  A,y(0)  by  A,y .  The  description  of  interest  to  us  is  then 


z(t )  =  D(e)Az(t )  (19a) 

or 


( ^(0  ^ 

(  4, 

0 

0  ^ 

(  An 

A12 

Aim  ^ 

( \ 

*■2  (*) 

— 

0 

eh2 

0 

A2i 

a22 

A2m 

^2(0 

VM*)y 

l  o 

0 

^  Am  i 

Am2 

a:  ) 

* - - 

The  objective  of  this  subsection  is  to  present  the  tests  of  MSSNS  and  MSST  for  the  reduced 
explicit  form  above.  It  is  easy  in  principle,  though  notationally  a  little  cumbersome,  to  apply 
the  operations  specified  in  Section  2.1  to  the  system  (19)  and  thereby  deduce  the  tests.  We 
shall  simply  state  the  results  obtained  in  [14].  One  might  anticipate  these  results  in  view  of 
the  similarity  of  (19)  to  (7);  the  conditions  for  MSSNS  (respectively  MSST)  are  those  that  one 
would  naturally  impose  in  order  to  get  an  eigenvalue  decomposition  (respectively  time-scale 
decomposition)  of  the  type  in  Example  1. 


To  begin,  suppose  the  x  kL  matrix  An  is  nonsingular.  Let  A2  denote  the  matrix 
obtained  as  the  Schur  complement  of  An  in  A(=  A(0)),  and  let  A22  denote  its  leading  k2  x  k2 
submatrix,  so 


Ao  — 


/  A22  . . .  A2m  ^ 


V  Afjx  2  ...  Am  m  J 


\ 

(  A21  } 

J 

^  Am  1  J 

Au*  (  Aj2  .  .  .  Aim  ) 


(20a) 


and 


A22  —  A22  —  A2i  A,  }  A 


lu 


12 


(206) 


Similarly,  if  A22  is  nonsingular,  let  A3  denote  the  Schur  complement  of  A22  in  A2,  and  let  A33 
denote  its  leading  k3  x  k3  submatrix.  The  pattern  in  this  construction  is  now  evident.  If  at 
some  stage  j  we  have  k}-  =  0,  then  the  corresponding  stage  is  skipped;  we  then  relabel  Ay  as 
Ay+i,  denote  its  leading  fcy+1  x  fcJ+1  submatrix  by  A}+U)+l,  and  proceed.  Note  that  Aj  is  to 
be  taken  as  A  and  An  as  An.  The  result  now  is  the  following: 

The  matrix  D(e)A  in  (19),  and  therefore  A(e)  as  well,  satisfies  MSSNS  (respectively 
MSST)  if  and  only  if  every  Ayy  for  which  fcy  0,  j  =  l,...,m,  is  nonsingular  (respec¬ 
tively  Hurwitz). 


(Note  that,  for  this  paper,  a  Hurwitz  matrix  is  one  whose  eigenvalues  are  in  the  open  left  half 
plane.)  Since  A  is  nonsingular  —  a  consequence  of  our  assumption  that  A(e)  is  nonsingular 
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away  from  0  —  the  nonsingularity  of  Amm  is  guaranteed  if  those  preceding  AJ}-  for  which  k}-  ^  0 
are  all  nonsingular.  The  index  m  here  can  be  seen  to  be  the  same  one  that  appears  in  Section 
2.1. 

It  was  mentioned  in  Section  2.2  that  P(e)  and  Q(e )  in  the  Smith  decomposition  (14)  are 
not  unique.  A  consequence  of  this  is  that  the  matrix  A  =  Q(0)P(0)  is  not  uniquely  determined 
by  the  decomposition.  However,  the  results  here  do  not  depend  on  which  decomposition  is 
chosen.  In  particular,  [14]  shows  that  the  above  Schur  complements  obtained  from  different 
decompositions  are  related  by  similarity  transformations. 
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3.  MSSNS/MSST,  EIGENVALUES  AND  INVARIANT  FACTORS 

3.1  MSSNS  and  the  Orders  of  Eigenvalues  and  Invariant  Factors 

This  subsection  will  firstly  establish  the  following  characterization  of  MSSNS: 

Result  1:  A(e)  has  MSSNS  if  and  only  if  the  orders  of  its  eigenvalues  equal  the  orders 
of  its  invariant  factors,  i.e.  if  and  only  if  it  has  precisely  kj  eigenvalues  of  order  j  —  1  for 
3  = 

An  analogous  result  is  presented  as  the  central  theorem  of  [5],  but  mapping  the  (high-gain 
feedback  root  locus)  formulation  and  proof  there  to  our  setting  here  would  lead  to  a  considerably 
more  cumbersome  proof  than  the  direct  one  below. 

In  the  process  of  proving  Result  1,  we  shall  also  be  demonstrating  the  earlier  claim  that 
under  MSSNS  the  eigenvalues  of  A(e)  fall  into  separate  groups,  with  fcy  eigenvalues  of  order  j—  1 
in  the  j'-th  group,  j  =  1,. . .  ,m,  and  that  the  eigenvalues  in  these  groups  can  be  approximated 
via  the  eigenvalues  of  reduced-order  matrices.  The  precise  result  will  be  stated  after  proving 
Result  1. 

Proof:  Since  the  unimodular  similarity  transformation  P(e)  described  in  Section  2.2  pre¬ 

serves  MSSNS/MSST,  as  well  as  eigenvalues  and  invariant  factors,  it  suffices  to  consider  a 
matrix  A(e)  that  is  already  in  explicit  form: 


A(e) 

D(e)A{e) 

(21a) 

(Ik, 

0 

0  A 

/  AnW 

Ai2fe) 

. .  Aim  (e)  ^ 

0 

0 

A21  (e) 

A22  (c) 

^2m  (c) 

(216) 

V  o 

0 

U-i(<o 

Am2(e) 

•  •  Amm  (e)  j 

with  A(e)  unimodular. 

To  show  the  “only  if”  part  first,  assume  that  A(e)  has  MSSNS.  Following  arguments  in 


[12],  let 

dy(e,  A)  =  det[e  1'A(e)-AJ] 

(22) 

and 

/y(e,A)  =  detF(e,A) 

(23a) 

where 

0  0  A 

0 

P~2Ik,  ■ 

F(e,X)  = 

l  0 

.  0 
'  ^  Ik j  +  kj+ 1  + ...+  km  J 

>"c 

■< 

\ 

—4 

1 

1 

(236) 
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so  that 


j'-i 

fi  (e> A)  =  «*  di  (<+  A) ,  v  =  J2  U  -  0  ki  (24) 

t  =  l 

For  e  ^  0,  dy(e,  A)  and  /y(e,  A),  regarded  as  polynomials  in  A,  will  have  the  same  roots.  Also, 
/y(e,  A)  is  a  continuous  function  of  e  at  e  =  0,  with 


/A  n 
(  A21 


f j  (0,  A)  =  det 


A12 

A22 

Ay  2 
0 


Aiy 

•  Ajm  ^ 

A2y 

A2  m 

*>■< 

..  1  , 
*v» 

Aym 

+  J  +  ...+  km  J 


(25) 


The  A^y  here  are  the  same  as  in  (19),  i.e.  submatrices  of  the  reduced  explicit  form.  The  latter 
determinant  is  easily  evaluated  by  iteratively  using  the  fact  that 


det  ( A2!  A22  )  =  det  Al  1  det  ^22  (26) 

if  An  is  nonsingular.  The  result  of  the  evaluation,  invoking  the  characterization  of  MSSNS  in 
Section  2.3  to  ensure  that  the  Schur  complements  A,-,-  below  (which  are  the  same  as  those  in 
Section  2.3)  are  nonsingular  for  *  =  1, . . .  ,j  -  1,  is 


J  1  m 

/y(0, A)  =  (H  det  A,-,)(det[Ayy  -  A i])A“,  tn  =  £>  (27) 

*=i  j+i 

Thus  /y(0,A)  has  k}-  roots  that  are  identical  to  the  kj  eigenvalues  of  Ayy,  which  we  shall  denote 
by  A ,(Ayy),  t'  =  1 , . . . ,  kj . 

It  follows  that  /y(e,  A)  and  hence  dy(e,  A)  have  kj  roots  that,  after  appropriate  matching, 
can  be  made  arbitrarily  close  to  these  respective  eigenvalues  by  choosing  e  small  enough.  Now 
(22)  shows  that  the  roots  of  dy(e,  A)  are  the  eigenvalues  of  A(e).  Hence  there  are  kj 

eigenvalues  of  A(e),  which  we  shall  denote  by 

1-1 

^»j+»>  si  ~  kg ,  t  =  1 , . . . ,  fcy  (28) 

9=1 


such  that  A  ,.+i/e3~1  can  be  brought  arbitrarily  close  to  A;(Ayy)  by  choosing  e  small  enough  (we 
are  assuming  for  notational  simplicity  that  the  indexing  follows  the  matching).  The  conclusion, 
since  all  the  eigenvalues  of  Ayy  are  nonzero  by  the  assumption  of  MSSNS,  is  that  A(e)  has 
precisely  fcy  eigenvalues  of  order  j  —  1 . 

For  the  “iP  part,  suppose  that  Ayy  is  singular,  but  A,»  for  t  =  1, ... ,  j  —  1  are  nonsingular. 
Then  /y(0,  A)  in  (27)  has  at  least  u>  +  1  roots  at  0,  so  A(e)  has  at  least  w  +  1  eigenvalues 

that  can  be  made  arbitrarily  close  to  0  by  choosing  e  small  enough.  It  follows  that  A(e)  has  at 
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least  w  +  1  eigenvalues  of  order  higher  than  j  —  1.  However,  it  has  only  w  invariant  factors  of 
order  higher  than  j  —  1.  •• 

This  result  is  of  intrinsic  interest  as  a  characterization  of  MSSNS.  It  also  leads  naturally  to 
the  results  of  Section  4  and  provides  approaches  to  the  problem  of  inducing  MSSNS  by  amplitude 
scaling  of  state  variables  via  (nonunimodular)  e-dependent  diagonal  similarity  transformations. 
The  proof  of  Result  1  also  establishes: 

Result  2:  Under  MSSNS,  the  eigenvalues  of  A(e)  converge  asymptotically  to  the  eigen¬ 
values  of  e,_1  Ajj,  j  =  1, . . .  ,m. 

3.2  Relating  MSST  to  MSSNS  and  the  Eigenvalues  of  A(e) 

It  is  of  interest  to  see  precisely  what  sorts  of  stability  conditions  on  the  original  system 
(l)  will  make  MSST  a  consequence  of  MSSNS.  The  new  result  to  be  proved  in  this  subsection 
is  the  following: 

Result  3:  A(e)  has  MSST  if  and  only  if  (i)  it  has  MSSNS,  and  (ii)  it  is  Hurwitz  for 
e  S  (0,  e0] ,  and  (iii)  the  order  of  the  real  part  of  every  eigenvalue  is  not  greater  than  the 
order  of  its  imaginary  part. 

The  third  condition  requires  that  the  order  of  the  damping  be  at  least  as  significant  as  the 
order  of  the  oscillation  frequency. 

Some  elementary  examples  may  help  to  illustrate  the  result  before  we  present  the  proof. 
The  first  condition  rules  out  the  matrix 

(o'  -e)  l29> 

from  having  MSST,  even  though  it  satisfies  the  other  two  conditions,  but 

(o'  -e)  (3°) 

satisfies  all  three  conditions  and  has  MSST.  The  third  condition  rules  out 

{-1  -e)  <3I> 

Proof:  To  prove  the  “only  iP  part,  note  from  Section  2.3  that  MSST  implies  that  all  the 

Ajj,  j  =  1  ,...,m,  are  Hurwitz.  They  are  hence  a  fortiori  nonsingular,  so  A(e)  has  MSSNS 
and  (i)  is  established.  It  then  follows  from  Result  2  that  the  eigenvalues  of  A(e)  converge 
asymptotically  to  those  of  so  (ii)  and  (iii)  follow. 
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For  the  “if”  part,  it  follows  from  (i)  and  Result  2  that  the  eigenvalues  A,  +t- ,  i  =  1 , . . . ,  Jby , 
j  —  1  of  A(e)  converge  asymptotically  to  those  of  e,-1Ayy.  From  the  nonsingularity  of 

the  Schur  complements  Ayy,  we  see  that 

lim(A,  i+i/e»'-1)  =  A, (A,-,-)  yt  0  (32a) 

It  follows  from  (iii)  that 

limRe  (A,i+t-/c,_x)  =  Re  (A<(Ayy))  yt  0  (326) 

Now  (ii)  implies  that  the  limit  above  is  negative,  and  it  follows  that  the  Schur  complements  are 
Hurwitz.  The  results  of  Section  2.3  then  show  that  A(e)  has  MSST.  •• 

The  assumption  of  MSST  has  dominated  studies  of  time-scale  decompositions.  However, 
systems  such  as  (31)  that  satisfy  conditions  (i)  and  (ii)  of  Result  3,  but  not  condition  (iii),  have 
been  considered  in  [14],  [15],  where  it  is  shown  that  so-called  extended  time-scale  decompositions 
can  be  obtained  under  appropriate  conditions.  The  basic  idea  here  is  to  retain  some  critical 
e-dependent  terms  in  the  scale-based  decomposition,  rather  than  simplifying  all  the  way  to  the 
constant  matrices  Ayy. 
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4.  MSSNS,  PRINCIPAL  MINORS  AND  AMPLITUDE  SCALING 


4.1  MSSNS  and  Principal  Minors 

Denote  the  eigenvalues  of  A(e)  by  A,-,  t  =  1 ,n,  and  its  j  xj  principal  minors  by  , 
i  =  1, . . . ,  ^  jQ  ,  j  =  1, . . . ,  n.  It  is  then  well  known  that 

det[A/  -  A(c)]  =  An  -  A-1  (£  A,)  +  A”-2 (£  \{\k)  -  . . .  +  (-1)"^  •  •  •  An)  (33a) 

iyik 

=  Xn  -  +  (336) 

Evidently  the  characteristic  polynomial  and  hence  the  eigenvalues  of  A(c)  are  determined  by 
its  principal  minors.  In  addition,  the  orders  of  the  coefficients  of  the  characteristic 

polynomial  suffice  to  determine  the  (fractional  or  integer)  orders  of  the  eigenvalues,  via  the 
classical  “Newton  polygon”  construction  reviewed  below,  see  also  [5],  [19]. 

The  invariant  factors  of  A(e)  are  determined  by  the  gcd’s  of  all  its  minors  of  each  di¬ 
mension.  However,  Result  1  and  (33)  lead  one  to  conclude  that,  under  MSSNS,  the  orders  of 
the  suffice  to  compute  the  invariant  factors.  We  show  in  this  subsection  that  under 

MSSNS  the  gcd’s  of  the  principal  minors  of  each  dimension  suffice  to  determine  the  invariant 
factors.  In  the  process,  we  clarify  the  role  that  MSSNS  plays  in  this  determination,  obtain  an 
alternate  characterization  of  MSSNS,  and  develop  insights  that  will  be  useful  in  the  amplitude 
scaling  problem  considered  in  the  next  subsection.  Again,  there  are  connections  with  results  in 
[5],  and  the  reader  may  wish  to  explore  these  (see  especially  Appendix  C  in  [5]). 

In  order  to  state  our  results,  some  additional  notation  is  needed: 

a{  will  denote  the  order  of  the  t-th  invariant  factor,  so  a{  —  j  —  1  if  j  is  the  smallest  integer 
for  which  i  <  ki  +  . . .  +  k}-,  with  the  kj  defined  as  in  Section  2.2; 

bi  will  denote  the  (possibly  fractional)  order  of  the  t-th  eigenvalue,  with  the  eigenvalues 
assumed  to  be  numbered  such  that  bi  <  ...  <bn; 

Pj  will  denote  the  order  of  the  gcd  of  all  j  xj  principal  minors  M of  A(e);  and 

r j  will  denote  the  order  of  the  sum  of  these  principal  minors. 

As  mentioned  above,  the  standard  Newton  polygon  construction,  [5],  [19],  can  be  applied 
to  the  relationships  embodied  in  (33a, b)  to  compute  the  b{  from  the  ry .  The  basis  for  this  lies 
in  the  following  lemma. 

Lemma  1:  Given  a  set  of  real  numbers  x,,  j  =  l,...,n,  there  is  a  unique  set  of  real 
numbers  y, ,  *  =  1, . . .  ,n  such  that  the  following  hold: 

yi  <  ••  •  <  y„  (LI) 
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y<  -  xi  when  yj  ^  y}+ 1 or  i  -  n  (i3) 

»=i 

The  numbers  t/<  may  be  obtained  as  the  slopes  of  the  segments  of  the  “lower  hull”  in  the 

(Newton  diagram)  plot  of  Xy  versus  j  (with  the  origin  included,  i.e.  x0  =  0),  see  Fig.  1. 

The  values  of  j  that  appear  in  (L3),  i.e.  values  where  t/y  ^  yy+i  or  j  =  n,  will  be  termed  corner 
points. 

Though  this  lemma  is  not  explicitly  articulated  in  the  usual  approaches  to  the  Newton 
polygon,  it  is  easily  seen  to  underlie  them.  Both  for  this  reason  and  because  a  proof  follows 
easily  from  the  construction  described  in  the  lemma  and  Fig.l,  we  omit  the  proof  here;  see  [14] 
for  details.  The  lemma  is  of  potential  value  in  other  settings  as  well,  which  is  why  we  have 
isolated  it. 

What  makes  Lemma  1  applicable  to  computing  the  6f  from  the  r3-  is  the  fact  that  these 
two  sets  of  numbers  satisfy  the  above  inequalities,  with  the  substitutions  xy  «—  ry  and  y,-  <—  6,-. 
To  see  this,  note  first  that  (Ll)  is  trivially  satisfied,  by  definition.  Also,  (33)  shows  that  the 
sum  of  the  j  x  j  principal  minors  equals  the  sum  of  all  j'-fold  products  of  the  eigenvalues,  so 
that  (L2)  is  satisfied.  Equality  fails  to  hold  in  (L2)  precisely  when  there  are  cancellations  of  the 
lowest  order  terms  when  summing  the  j-fold  products  of  eigenvalues.  For  example,  if  n  =  2, 
Ai  =  -1  and  A2  =  1  +  e,  then  =  0  (and  i>2  =  0)  but  tq  =  1.  However,  if  b}-  ^  b]  +  1  then  there 
will  be  only  one  j-fold  product  of  eigenvalues  that  has  order  bi  +  . . .  +  bj,  namely  the  product 
Ai  •  •  •  Ay,  and  all  other  j-fold  products  will  have  higher  order,  so  equality  holds  in  (L2).  Also,  if 
j  =  n  then  there  is  only  one  term  in  the  sum,  namely  Ai  •  •  •  A„,  so  equality  holds  again.  Thus 

(L3)  is  established.  The  conclusion  is  that  the  b{  are  given  by  the  slopes  of  the  segments  of  the 

lower  hull  in  the  Newton  diagram  for  the  points  ry . 

Some  further  connections  between  the  numbers  defined  above  should  be  noted.  In  general 
it  will  be  the  case  that 

Pi  <  ri  (34) 

Equality  fails  to  hold  precisely  when  the  leading  (lowest-order)  terms  cancel  out  upon  summing 
the  j  x  j  principal  minors.  When  equality  holds  in  (34),  we  shall  say  that  the  no-cancellation 
condition  holds.  It  is  also  evident  from  the  definition  of  the  a,  that 

a i  <  ...  <an  (35) 


Furthermore,  the  definition  of  invariant  factors  in  Section  2.2  shows  that  ai  +  . . .  +  ay  is  the 


order  of  the  gcd  of  all  j  x  j  minors  of  A(e),  so  in  particular 

j 

Y,a<^Pi  (36) 

t=  1 

What  (35)  and  (36)  demonstrate  is  that,  with  the  substitutions  x}-  *—  py  and  yt-  *—  ay,  the 
integers  py  and  a,-  satisfy  (Ll)  and  (L2).  If  (L3)  was  also  satisfied,  then  Lemma  1  would  imply 
that  the  a,-  can  be  uniquely  determined  from  the  py ,  as  the  slopes  of  the  segments  of  the  lower 
hull  in  the  Newton  diagram  for  the  points  py. 

The  next  result  shows  that  MSSNS  guarantees  (L3)  with  the  above  substitutions,  so  under 
MSSNS  the  invariant  factors  can  indeed  be  determined  from  the  gcd’s  of  the  principal  minors  of 
each  dimension,  using  the  Newton  polygon  construction  associated  with  Lemma  1.  The  result 
below  also  shows  that  under  MSSNS  (34)  holds  with  equality  (i.e.  the  no-cancellation  condition 
holds)  at  the  corner  points.  Conversely,  if  (L3)  holds  with  the  preceding  substitutions,  and  if 
the  no-cancellation  condition  holds  at  the  corner  points,  then  A(e)  satisfies  MSSNS.  When  A(e) 
is  in  explicit  form,  then  the  no-cancellation  condition  need  not  be  checked:  A(e)  in  explicit  form 
has  MSSNS  if  and  only  if  (L3)  holds  with  the  preceding  substitutions. 

Result  4:  A(e)  has  MSSNS  if  and  only  if 

i 

X>f  =Py  =  ry  when  Oy^ay+1  (37) 

»'=  1 

Furthermore,  when  A(e)  is  in  explicit  form,  only  the  first  equality  in  (37)  needs  to  be 

checked,  because  it  implies  the  second. 

Proof:  To  prove  the  “only  if”  part  first,  assume  A(e)  satisfies  MSSNS.  Then  from  Result 

1  we  have  that  b{  =  a,-,  t  =  1, . . .  ,n.  Now  suppose  ay  yt  ay+1.  Then  6y  yt  bJ+1,  so  by  our  earlier 
arguments  bi  +  ...-(-  b}-  =  ry.  Hence  ai  +  . . .  +  ay  =  ry .  Combining  this  with  the  inequalities  in 
(34)  and  (36)  demonstrates  the  result. 

For  the  “if”  part,  note  from  (34)  that  the  lower  hull  for  the  points  ry  in  the  Newton 
polygon  construction  cannot  lie  below  that  for  the  py.  Furthermore,  since  (37)  shows  that  the 
corner  points  for  the  ry  hull  coincide  with  those  for  the  py  hull,  the  two  hulls  must  be  the 
same.  Hence  the  numbers  6t-  obtained  by  applying  Lemma  1  to  the  ry  must  be  the  same  as  the 
numbers  a,-  obtained  by  applying  Lemma  1  to  the  py.  In  other  words,  6,-  =  a<5  so  by  Result  1 
the  matrix  A(e)  has  MSSNS. 

Suppose  now  that  A(e)  is  in  the  explicit  form  (16)  and  that  the  first  equality  in  (37)  holds, 
so  ay  yt  aJ  +  1  and  ax  +...-(-  ay  =  py .  Now  examination  of  the  explicit  form  shows  that  there  is 
only  one  principal  minor  of  order  py ,  namely  the  leading  principal  minor  of  appropriate  size  (its 
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size  is  actually  ki  + . . .  +  fc,  where  q  is  the  smallest  integer  for  which  k2  +  2 fc3  + . . .  +  qkv+  x  >  py  ) . 
Hence  r}-  =  py,  i.e.  the  second  equality  in  (37)  holds  and  the  no-cancellation  condition  is  satisfied 
at  the  corner  points.  •• 

The  need  to  check  the  second  equality  in  (37)  when  A(e)  is  not  in  explicit  form  can  be 
illustrated  by  the  case  of  the  matrix 

(1  -r+\0  <38> 

which  has  px  =  0,  p2  =  2  and  ax  =  0,  a2  =  2  so  (L1)-(L3)  are  satisfied  with  the  substitutions 

above,  but  rx  =  2,  r2  —  2  and  6X  =  1,  b2  =  1,  so  indeed  the  matrix  does  not  have  MSSNS.  For 

an  example  to  illustrate  how  we  might  use  all  the  above  results,  consider  the  following. 

Example  2  Let 


It  is  easily  seen  that  this  is  not  in  explicit  form.  To  determine  the  a{ ,  we  use  the  definition  of 
invariant  factors  in  Section  2.2,  which  involves  examining  all  minors  of  each  dimension.  A  quick 
inspection  shows  that  there  are  lxl  and  2x2  minors  that  are  of  order  0,  so  ax  —  0,  a2  —  0. 
There  is  a  3  x  3  minor  (which  happens  to  be  the  leading  principal  minor)  of  order  3,  and  none 
of  order  2,  hence  a3  =  3.  Finally,  the  determinant  is  of  order  10,  so  a4  =  7.  Examining  the 
principal  minors,  we  see  that  px  =  1,  p2  =  0,  p3  =  3,  p4  =  10.  Hence  the  a{  and  py  satisfy 
(L1)-(L3),  with  the  former  being  the  slopes  of  the  lower  hull  of  the  latter  set  of  points  in  the 
Newton  diagram.  Now  if  the  no-cancellation  conditions  hold  at  the  corner  points,  i.e.  if  r2  =  p2 
and  r3  =  p3,  then  we  will  be  able  to  conclude  from  Result  4  that  A(e)  does  indeed  satisfy 
MSSNS.  Since  there  is  only  one  2x2  principal  minor  of  order  0  and  only  one  3x3  principal 
minor  of  order  3,  these  conditions  are  indeed  satisfied.  With  the  assurance  that  A(e)  satisfies 
MSSNS,  one  can  now  proceed  with  the  additional  work  required  to  transform  it  to  explicit  form 
and  carry  out  scale-based  decompositions  of  it.  •• 

4.2  Amplitude  Scaling  to  Induce  MSSNS 

While  unimodular  similarity  transformations  do  not  affect  MSSNS  (or  MSST),  nonuni- 
modular  similarity  tranformations  may  do  so,  by  modifying  the  invariant  factors  of  the  matrix 
(the  eigenvalues  are  of  course  preserved).  In  particular,  such  transformations  may  be  used  to 
induce  MSSNS  in  a  matrix  that  does  not  satisfy  it,  after  which  scale-based  decompositions  may 
be  carried  out  as  earlier.  Our  aim  in  this  section  is  only  to  illustrate  what  is  possible. 


We  restrict  ourselves  to  nonunimodular  transformations  that  are  diagonal,  which  corre¬ 
spond  to  e-dependent  amplitude  scaling  of  the  individual  state  variables.  Such  transformations, 
in  addition  to  preserving  eigenvalues,  also  preserve  all  principal  minors  of  the  matrix  they  act 
on.  The  results  of  the  previous  subsection  then  show  that  a  necessary  condition  for  such  a 
transformation  to  induce  MSSNS  is  that  the  no-cancellation  condition  p3  =  r,  be  satisfied  at 
the  corner  points  of  the  Newton  diagram  for  either  of  these  sets  of  indices. 

Our  amplitude  scaling  results  are  drawn  from  the  thesis  [14],  which  includes  results  on  a 
systematic  amplitude  scaling  procedure  for  matrices  A(e)  that  satisfy  certain  conditions.  Since 
these  conditions  are  not  only  rather  strong  but  also  hard  to  test  for,  we  do  not  attempt  to  do 
more  than  illustrate  a  simple  case  of  the  procedure  here.  We  begin  with  a  canonical  example 
that  serves  to  illustrate  the  basic  idea  behind  the  procedure. 


Example  3  Let 


A(e)  = 


( 0 

0  . 

0 

\ 

0 

0 

€°a  . 

0 

0 

0 

0  . 

.  ean~ 

0 

0  . 

0 

/ 

(40) 


The  nonnegative  integers  a,-  can  be  seen  to  be  the  orders  of  the  invariant  factors  of  A(e),  though 
not  necessarily  in  ascending  order  now.  The  j  x  j  principal  minors  for  j  <  n  and  hence  the 
corresponding  py,  r}-  are  all  0,  while  a1  +  . . .  +  an  =  pn  =  rn.  By  applying  the  Newton  polygon 
construction  to  the  ry  or  by  directly  evaluating  the  characteristic  polynomial  of  A(e),  it  is  easily 
seen  that  the  eigenvalues  all  have  the  same  order,  namely  b  =  rn/n.  Thus  A(e)  will  have  MSSNS 
if  and  only  if  all  the  a<  are  equal  (and  equal  to  b). 

Now  defining  the  similarity  transformation 


S (e)  =  diagonal  [e1" 1 , . . . ,  ew  “_1 , 1] 


(41a) 


with 

=  «><+i  +  b-a{, 

i  = 

i,... 

,,n-  1, 

Wn  - 

=  0 

(416) 

it  is  easily  verified  that 

(0 

eb 

0  ... 

0 

0 

e6  ... 

0 

S(e)A{e)S~'(e)  = 

0 

0 

0  ... 

(42) 

U" 

0 

0  ... 

0  J 

so  that  this  transformed  matrix  has  MSSNS.  •• 


The  following  procedure  is  suggested  by  examples  such  as  the  above,  and  can  be  guar¬ 
anteed  to  work  under  certain  strong  conditions,  [14].  The  first  step  is  to  transform  A(e)  to  its 
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explicit  form.  We  assume  from  now  on  that  this  has  been  done,  and  use  A(e)  to  denote  the 
matrix  in  explicit  form.  In  the  second  step,  we  identify  what  may  be  termed  a  skeleton  in  A(e). 
A  skeleton  consists  of  n  elements,  precisely  one  from  each  row  and  column  of  the  matrix,  with 
the  order  of  the  element  in  the  t-th  row  being  a,- .  There  has  to  be  at  least  one  skeleton  in  A(e) 
because  of  our  standing  assumption  of  nonsingularity  away  from  0. 


In  the  third  step,  we  similarity  transform  A(e)  with  a  permutation  matrix  that  brings 
the  elements  of  the  skeleton  to  the  locations  of  the  l’s  in  the  block  diagonal  canonical  circulant 
matrix,  whose  diagonal  blocks  take  the  form 


/  0  1  0  ...  (A 
0  0  1  ...  0 


0  0  0  ...  1 
ll  0  0  ...  0  J 


(43) 


Though  [14]  considers  the  case  of  multiple  blocks,  we  restrict  ourselves  here  to  the  case  of  a 
single  block,  i.e.  to  the  case  where  the  elements  of  the  skeleton,  after  transformation,  lie  at 
the  locations  of  the  n  x  n  circulant  matrix  (43).  Let  a,-  now  denote  the  order  of  the  skeleton 
element  in  the  t-th  row  of  the  transformed  matrix. 

The  final  step  of  the  procedure  is  to  transform  the  matrix  with  the  similarity  transfor¬ 
mation 

5(e)  =  diagonal  [e*"1,...  ,e'"n_l,l]  (44a) 


with 


Wi  =  tvi+ 1  +  i>i  -  a,- ,  t  =  1, . . . ,  n  -  1,  wn  —  0 


(446) 


where  the  b{  are  the  orders  of  the  eigenvalues.  Under  conditions  described  in  [14] ,  the  resulting 
matrix  satisfies  MSSNS. 

We  illustrate  this  procedure  with  an  example. 


Example  4 
in  brackets: 


Consider  the  explicit  form  matrix  below,  with  the  skeleton  elements  enclosed 


A(e)  = 


fe3 

€4 

e5 

un 

[e] 

e3 

e3 

€ 

e3 

£2 

e7 

\e6 

£8 

[£«] 

<?) 

(45) 


Since  the  matrix  is  in  explicit  form,  it  is  easy  to  check  for  MSSNS.  The  matrix  An  referred  to 
in  Section  2.3  is  simply  the  (1,1)  entry  of  the  A(e)  here  evaluated  at  e  =  0,  and  is  0,  so  MSSNS 
does  not  hold.  This  fact  can  also  be  seen  after  determining  that  ax  =  0,  a2  =  1,  a3  =  1,  a4  =  6 
and  6i  =  62  =  63  =  64  =  2. 
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The  similarity  transformation  that  brings  the  skeleton  elements  to  the  canonical  positions 
mentioned  above  simply  involves  interchanging  the  first  and  third  rows  and  then  the  first  and 
third  columns.  The  result  is  the  matrix 


AiM 


(  e2  [«]  e3  eM 
e3  e3  [e]  e 
e4  [i] 

W)  *  36  «T; 


(46) 


Now  using  (44)  we  find  that  twi  =  4,  w2  =  3,  w3  —  2,  tw4  =  0.  Transforming  Ai(e)  with  the 
resulting  S(e)  gives 


A2  ( e )  =  5  (e)  Ax  (e)S  1  (e)  = 


/  e2 

e2 

e5 

e11^ 

e2 

e3 

e2 

e4 

e3 

e3 

e2 

e2 

V2 

e3 

e4 

<?  ) 

(47) 


which  is  easily  seen  to  have  MSSNS  (it  is  in  explicit  form,  with  invariant  factor  orders  all  equal 
to  2).  •• 


The  amplitude  scaling  procedure  illustrated  above  is  further  developed  in  [14].  It  has 
been  used  to  motivate  the  amplitude  scaling  carried  out  in  [20],  and  to  treat  some  further 
generalizations. 
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5.  CONCLUSION 


The  algebraic  approach  to  the  analysis  of  singularly  perturbed  systems  of  the  form  (1), 
as  described  here  and  in  [14],  [15],  has  added  a  useful  dimension  to  what  continues  to  be  an 
active  area  of  study.  It  has  provided  novel  ways  to  understand  and  interpret  previous  results, 
and  has  yielded  new  results  and  new  tools  for  analyzing  singularly  perturbed  systems.  In 
particular,  the  MSSNS  and  MSST  conditions  that  have  been  crucial  to  earlier  work  in  singular 
perturbations  have  been  examined  and  illuminated  from  the  algebraic  viewpoint  in  this  paper. 
The  connections  with  (frequency-domain)  techniques  and  results  developed  in  the  study  of 
asymptotic  root  loci  are  also  much  more  evident  in  our  setting. 

Fresh  questions  have  arisen  naturally  in  the  context  of  this  algebraic  approach.  For 
example,  it  is  of  interest  to  further  understand  conditions  and  procedures  for  the  extended 
time-scale  decompositions  mentioned  briefly  at  the  end  of  Section  3.  Also,  Section  4  has  only 
touched  on  the  task  of  inducing  MSSNS  by  nonunimodular  similarity  tranformations,  and  has 
not  discussed  in  any  detail  the  problem  of  relating  properties  of  the  transformed  system  back 
to  the  original  system.  Another  worthwhile  task  is  to  understand  connections  with  the  Jordan 
structure  of  A(e),  perhaps  making  deeper  contact  with  the  results  of  [10],  [21]  in  the  process. 
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Fig.  1.  Newton  polygon  construction. 
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